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Quarkonium production is a sensitive probe of the dynamics of hard scattering, which can 'measure' 
the environment of the heavy quark pair after it is created in a hard process. Quarkonium hadropro- 
duction data indicates that the produced pair experiences a secondary, hard interaction and then 
ceases to interact at a stage when the pair is still compact (compared to the size of the quarkonium 
wave function). These features differ from scenarios mostly discussed so far. An approach which 
relies on an early, perturbative rescattering of the pair with a comoving color field can explain many 
observed features, including the polarization, the X1/X2 ratio and the nuclear target ^4-dependence, 
which are difficult to understand otherwise. 
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I. INTRODUCTION 

The production of charm and bottom quarks provides 
a rich testing ground for the dynamics of hard QCD pro- 
cesses. The detection of heavy quark(s) is an unambigu- 
ous signal of a hard underlying process at the quark mass 
scale niQ, with no need for jet reconstruction. The heavy 
quark momentum is (roughly) determined by that of the 
heavy hadron in which it is contained. The quarks are 
colored and reinteract with other colored partons in the 
final state, making heavy quark processes complementary 
to lepton pair (Drcll-Yan) production. Existing data on 
photo- and hadroproduction of open and hidden charm 
and beauty production, at small and large transverse mo- 
menta and on a variety of targets allows systematic stud- 
ies of the properties of hard processes. For a recent review 
of many aspects of open heavy quark production see Ref. 

In the standard QCD approach to hard inclusive pro- 
cesses the cross section is written as a product of uni- 
versal structure and fragmentation functions times a 
process-specific PQCD subprocess. This is justified by 
the QCD factorization theorem [gj, which applies when 
the sum over the final hadronic states is well approxi- 
mated by the partonic one. For heavy quarks this is typ- 
ically the case well above production threshold. In the 
total heavy quark cross section the average transverse 
momentum of a heavy quark is comparable to its mass, 
p± ~ toq. For large wiq many final states are avail- 
able, and the heavy quark hadronization is kinematically 
well separated from the spectator partons of the colliding 
hadrons. The total charm and beauty cross sections in 



both photo- and hadroproduction in fact agree well with 
PQCD expectations Q. 

A closer look at the production characteristics of charm 
quarks nevertheless reveals some surprising features. One 
would expect the heavy quarks to be produced nearly 
back-to-back in azimuthal angle, i.e., they should be 
coplanar with the projectile. As seen from Fig. 1, this is 
rather well satisfied for charm photoproduction, whereas 
for hadroproduction the measured azimuthal distribution 
is nearly flat. To fit the hadroproduction data requires a 
'primordial' transverse momentum of the partons in the 
colliding hadrons as large as (kj_) ~ 1 GeV 2 , of the same 
order as the scale m;? of the hard process. 

The different characteristics of photo- and hadropro- 
duction suggests that the reason for the acoplanarity lies 
in gluon radiation of the projectile parton. This radi- 
ation creates a 'hot' environment for the heavy quarks 
which may further increase their acoplanarity through 
secondary interactions. Such interactions do not affect 
the factorization theorem for observables (like the total 
cross section) which involve an inclusive sum over the 
heavy quark momenta. 

In quarkonium production the relative momentum of 
the heavy quarks is severely constrained - to a small frac- 
tion of m,Q for a non-relativistic bound state - hence the 
QCD factorization theorem does not apply. There is no a 
priori reason to believe that the production cross section 
can be expressed as a product of universal parton distri- 
butions and a PQCD subprocess. In particular, the hot 
environment observed in DD production will be experi- 
enced also by the heavy quarks that bind into quarkonia. 
This environment may well affect the quarkonium cross 
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Solid: QCD NLO 

Dashed: NLO+k T kick, <k 2 > = 0.5 GeV 2 
Dotted: NLO+k T kick, <k?> = l GeV 2 
x : E687 data (yN) 
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FIG. 1. The relative azimuthal angle distribution of DD 
pairs in photoproduction (upper figure) and hadroproduction 
(lower figure). Figures taken from Ref. [1]. 



section - the bound state can be broken up by consid- 
erably softer interactions than required to get the nearly 
flat azimuthal distribution of Fig. 1. 

Here I would like to focus on recent work done together 
with S. Peigne j|. Based on properties of the quarkonium 
data we construct a scenario for the production dynamics 
which is consistent with QCD. We also discuss an explicit 
model which appears to explain several puzzling features 
of the data. 



II. HINTS FROM THE QUARKONIUM DATA 

As noted above (Fig. 1), there is evidence from DD 
production that hadroproduced heavy quarks are created 
together with a strong color field. If this field is causing 
the 'anomalies' of the quarkonium cross section then 
one would expect no anomalies in quarkonium photopro- 
duction, in the fragmentation region of the photon (which 
should be free of color fields). As a matter of fact, the 



Color Singlet Model (CSM) provides a very good descrip- 
tion of inelastic J/ip photoproduction . 

The time scale of reinteraction effects between the 
heavy quarks and the surrounding color field in hadropro- 
duction can be surmised from the ratio a(ip')/adi r (J/ip) 
of the ip' and the directly (i.e., not through cascade de- 
cays) produced J/ip. The measured ratio is known to 
be independent of the nature of the projectile and target 
and of the kinematics. Its value is furthermore consistent 
with the square of the ratio of the wave functions at the 
origin P-pT| . This suggests that the reinteractions occur 
while the heavy quark pair is still in a compact config- 
uration, compared to the spatial size of the quarkonium 
wave function. 

A further hint concerning the quarkonium produc- 
tion dynamics is provided by the nuclear target A- 
dependence. The J/ip and ip' cross sections both scale 
as a A ^ cr N A a , with a = 0.92 ± 0.01 jnj. On the other 
hand, for open charm (D) production a — 1.02 ± 0.03 ± 
0.02 p^ |. This difference in A-dependence is quite re- 
markable. At high energies the charm pair remains in 
a compact ~ l/m c configuration while it is inside the 
nucleus. The uncertainty in its energy is thus so large 
that the pair couples to both open and hidden charm 
channels. In other words, the 'decision' as to whether 
a given quark pair will materialize as a J/ip, ip' or DD 
final state can be made only after the pair has left the 
nucleus. Apparently the only way of understanding their 
differing A-dependence is then that the formation of J/ip 
and ip' , but not of DD, requires a further interaction with 
comoving partons. If the comovers are 'swept away' by 
the nucleus a suppression of the quarkonium cross section 
will result @. 

The transverse separation of such comovers from the 
quark pair can be estimated from the effective 'absorp- 
tion' cross section needed to describe the observed nu- 
clear suppression. In a Glauber framework one finds 
a a b s — 7.3 ± 0.6 mb |14| which (using a = irr 2 ) corre- 
sponds to a radius of about 0.5 fm. This is considerably 
bigger than would be expected for a compact cc state 
(r = l/m c ~ 0.13 fm), although it is still smaller than 
the typical QCD scale 1/Aqcd — 1 fm. 

The anomalous effects in the cross section do not ap- 
pear to vanish quickly with the quark mass. The T(3S f ) 
cross section exceeds the CSM prediction by an order 
of magnitude, and the nuclear suppression persists (al- 
beit at a reduced level). Whatever causes the anomalies 
seems to be of 'leading twist', in the sense that its effect 
vanishes slowly with tuq . 

Similarly, the anomalous effects observed in the to- 
tal quarkonium cross section do not decrease (if any- 
thing, they increase) with the transverse momentum of 
the quarkonium. Thus, for p± 3> m c the CSM underesti- 
mates the measured J/ip and Xci cross sections by more 
than an order of magnitude, whereas the relative produc- 
tion rate of J / ip and ip' is close to expectations based on 
the wave function at the origin. Since the anomalous ef- 
fects are similar in the total and large p± cross sections it 



2 



is reasonable to assume that they have the same origin. 
The following discussion concerns the total cross section, 
but an analogous scenario is applicable to the large p± 
cross section. 



III. RESCATTERING FROM A COMOVING 
COLOR FIELD 

How can we understand and describe a reinteraction 
mechanism, qualitatively indicated by the data as dis- 
cussed above, in terms of QCD? The loop diagram shown 
in Fig. 2b may seem like a good candidate. The projectile 
gluon radiates a gluon which reinteracts with the quark 
pair. However, on closer inspection this diagram does 
not describe a two step process, with the reinteraction 
clearly separated in time and space from the production 
of the heavy quark pair. The spatial size of the loop is 
dictated by the hard scale to be l/m c . The reinterac- 
tion thus occurs simultaneously with the creation of the 
heavy quarks, and the diagram must be viewed as a ver- 
tex correction to the lowest order CSM diagram of Fig. 
2a. 



(2a) 




(2b) 
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FIG. 2. Basic processes for J/tp hadroproduction in the 
CSM (figures (2a) and (2b)) and in the present scenario (figure 
(2c)). 

A qualitatively different possibility is sketched in Fig. 
2c. Here the quark pair is first created in a color octet 
state, together with the color field radiated from the col- 
liding gluons. This system propagates for a time com- 



mensurate with the time-scale of emission of the radi- 
ated field (which is long compared to the time l/m c in 
which the quark pair is created) . During the propagation 
the emitted field interacts with itself and also with the 
quark pair. Most of these interactions are relatively soft 
'monopole' ones, where the quark pair acts as a pointlikc 
colored object. The internal quantum numbers of the 
pair (e.g., the quark spins) are thus unchanged. Even- 
tually a harder dipole interaction (gluon T in Fig. 2c) 
resolves the internal structure of the pair and changes it 
into a color singlet (in this process the quark spins may 
flip). If the color field has lost its phase coherence with 
the quark pair it will act like a classical source. 

The existence of a phase incoherent color field T as in 
Fig. 2c seems not to be excluded by any theoretical ar- 
guments. The qualitative comparison of QED and QCD 
production processes in Fig. 3 may help to explain how 
such a field could arise in hadroproduction, but still be 
absent in photoproduction. 

The upper half of Fig. 3 is a sketch of the QED pro- 
cess e + e~ — » /J + /i~ in its CM, before and just after the 
creation of the muon pair. The incoming electrons carry 
a radiation field, which is released when the electrons an- 
nihilate. The right- and leftmoving fields do not interact 
and thus simply pass each other, leaving the muon pair in 
an initially bare state, without a surrounding (comoving) 
field. 

In the QCD process gg — > QQ (lower half of Fig. 3) the 
colliding gluons again carry a radiation field. However, 
now the right- and leftmoving fields interact as they pass 
each other, creating field components at lower rapidities. 
Qualitatively, the situation seems analogous to the cre- 
ation of a rapidity plateau in hadron collisions. The final 
color field could thus be evenly distributed in rapidity, 
as sketched in Fig. 3. The external field T of Fig. 2c 
would then correspond to the central parts of this field, 
with low CM rapidity. 

In the subprocess 73 — > QQ of heavy quark photopro- 
duction, only the gluon carries a radiation field. After 
the creation of the heavy quarks this field then contin- 
ues without interactions and emerges at high rapidity in 
the final state. The bare heavy quarks are left without a 
comoving field (r = 0), much like in the QED example 
of Fig. 3. Hence one would expect the CSM to work in 
photoproduction, as observed experimentally. 



IV. AN EXPLICIT MODEL OF QUARKONIUM 
HADROPRODUCTION 

The scenario for quarkonium hadroproduction dis- 
cussed in the previous section is qualitative. Detailed 
predictions depend on the properties of the comoving 
field r. Here I would like to briefly describe the results 
of a model study, based on a number of simplifying as- 
sumptions ||. In particular, we assumed that T is weak 
enough so that only a single 'dipole' interaction needs to 
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FIG. 3. Upper half: Sketch of the e e~ — > p, fi~ pro- 
cess before (left) and just after (right) the creation of the 
muon pair. The motion of the radiation field (grey) is in- 
dicated by double arrows. Lower half: Analogous sketch of 
99 —* QQ- Interactions between the colliding colored radia- 
tion fields create field components which are comoving with 
the heavy quarks, i.e., which have low rapidity in the QQ rest 
frame. 



be considered, and that T is isotropically distributed in 
the rest frame of the QQ. 

The scattering amplitude consists of three parts, and 
is calculated from the diagram shown in Fig. 4. The 
first part is the gg — > QQ amplitude <E>, calculated using 
PQCD. The 'monopole' interactions which occur after 
the fusion process do not change the quantum numbers of 
the quark pair (except for its octet color component) and 
thus need not be explicitly taken into account. However, 
in the time interval between the creation of the pair and 
its reinteraction with T the spatial size of the quark pair 
increases by some factor p, which is a parameter of our 
model. 




FIG. 4. A perturbative interaction between the quark pair 
and a gluon from the color field F creates an overlap between 
the QQ wave function $ from the gg — > QQ process with 
the physical quarkonium wave function ^. There is a second 
diagram where the gluon attaches to the antiquark. 



The reinteraction of the quark pair with the external 
field r constitutes the second part of the amplitude of 
Fig. 4. Since this is a relatively hard 'dipole' interaction 
the rescattering amplitude TZ is again evaluated using 
PQCD. In the QQ rest frame the 3-momentum part of 
the transfer I dominates, <C \l\ <C TOq, since the 
heavy quarks are nearly on-shell both before and after 
the rescattering. 

The third part of the scattering amplitude is the 
quarkonium wave function If at the time of rescat- 
tering the quark pair is still compact compared to the 
size of physical quarkonium, the amplitude depends only 
on the quarkonium wave function at the origin (or on 
its derivative). For S-wave quarkonia the expansion fac- 
tor p is then irrelevant. For P-wave quarkonia, on the 
other hand, the overlap is oc p since their wave function 
vanishes at the origin. 

The full amplitude of Fig. 4 is thus expressed as a 
double convolution, 



M = ]T {LL Z ;SS Z \JJ Z ) ^ 



XX, aa 



d 3 p d 3 q 



X$g(p) n x - KaS {p,q)^> 



(2tt) 3 (2tt) 3 
''*(<?)* (1) 



Explicit expressions for the various parts of this ampli- 
tude may be found in || . It turns out that Eq. (Q) implies 
two distinct mechanisms for S'-wave quarkonium produc- 
tion, depending on which Lorentz components of the ex- 
ternal field T M are involved: 

(i) The creation of the QQ pair in an S = L = state, 
followed by a spin-flip interaction with a trans- 
versely polarized gluon from the color field T. 

(ii) The creation of the pair in an S = L = 1 state, fol- 
lowed by a spin-conserving interaction with a lon- 
gitudinally polarized gluon from the color field T. 

Case (i) leads to production of unpolarized J/i/j's, since 
information about the collision axis is not transmitted 
by the S = L = quark pair, and we assume T to be 
isotropic. On the other hand, mechanism (ii) gives rise 
to transversely polarized J/ip's. 

Experimentally [fl5| 17 both the J/tp and the if)' are 
unpolarized (at the moderate values of xf which are rel- 
evant for our scenario). We conclude that our model 
correctly describes the ip' polarization data provided the 
rescattering involves transverse gluons (which implies 
r° = 0). Conclusions for the J/tp are complicated by 
the fact that about 40% of the measured J/ ip cross sec- 
tion stems from Xci.2 — > J/tp + j radiative decays p8| , ^9| . 
This component has not been removed in the polarization 
measurements. Hence we need to take into account also 
the production and polarization of P-wave quarkonia. 

It turns out that in our model the P-wave cross sections 
have the same dependence on the external field T as the 
S'-wave ones. The cross section ratios are given in Table 
I in terms of 
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TABLE I. Relative 3 Pj cross sections and induced J/ip 
polarizations. 



3 ,/i£/mV /2.510-V (cc) 
r "5 P ^"rTJ -\6.510- 3 p 2 (66) l j 

which depends on the expansion paramter /? of the quark 
pair wave function. 

The total cross section ratio er(xi) : <j(x2) = 3:5 
of Table || is consistent with the experimental number 
0.6 ± 0.3. It should be kept in mind, however, that the 
Xc2 cross section calculated in the CSM is not far be- 
low the measured one. The CSM could well account 
for half of the measured <r(xc2)- Adding this contribu- 
tion to ours will correspondingly decrease the calculated 
o-(Xi) : fJ (X2) ratio. 

The measured c(Xc2) is somewhat larger than 
0dir{J/ip)- This is generically difficult to understand in 
an approach like ours, which couples a compact quark 
pair directly to the quarkonium wave function. As seen 
from Eq. (g), the P-wave cross section is suppressed due 
to the wave function vanishing at the origin. We need 
to have a considerable expansion of the quark pair wave 
function between its creation and reinteraction, p ~ 3, 
to describe the data. Such a large expansion makes the 
cc radius comparable to that of the J ftp, i.e., the quark 
pair is no longer 'compact'. 

In Table | we also show the induced J/ip polarizations 
from the radiative decays of the various P-wave states. 
The J/ip polarization is parametrized through the ratio 
A defined by 

= <r(J z = +1) - a(J z = 0) 

<j{Jz = +l) + a{J z = 0) [ 1 

The P-wave helicity states which are produced in our 
model induce a longitudinal (A < 0) polarization for the 
J /if). On the other hand, the CSM produces Xc2's with 
J z = ±2 ||, giving transversely polarized J/ip's. If the 
CSM contributes ~ 50% of the measured cr(xc2), the to- 
tal indirectly produced J/tp's turn out to be unpolarized. 
Since in our mechanism the directly produced J/ip com- 
ponent is also unpolarized, the experimentally observed 
non-polarization of the J / ip can thus be understood. 

V. CONCLUSIONS 

Quarkonium cross sections are sensitive to rescatter- 
ing between the produced heavy quarks and accompa- 



nying colored partons. The standard QCD factorization 
theorem which allows one to ignore such effects in hard 
inclusive processes does not apply to quarkonium pro- 
duction. The observed angular (de)correlations between 
hadroproduced DD pairs show that multiple scattering 
effects are strong. 

The quarkonium data as well as theoretical arguments 
suggest H a rescattering scenario where the relevant 
rescattering (i.e., the one which affects the cross section) 
is perturbatively hard and occurs soon after the heavy 
quarks are produced. The rescattering is associated with 
the gluon-gluon fusion process and thus depends on the 
colliding hadrons only via their universal gluon distribu- 
tions. There is no rescattering in the photon-gluon fusion 
process. 

Such a scenario is consistent with the phenomenolog- 
ical successes of the Color Evaporation Model (CEA4) 
p0| , pl| . It avoids a theoretical inconsistency of the 
CEM, namely that heavy quark spins would have to be 
flipped by soft interactions. It also avoids some phe- 
nomenological problems of the CEM, namely that ra- 
diative Xc decays are observed to contribute less in J/ip 
photoproduction than in hadroproduction p2| , and that 
quarkonium hadroproduction on nuclear targets is sup- 
pressed. Finally, the present scenario explains why the 
ratio <j(ip')/(Jdir(J/ip) is close to the ratio of the wave 
functions at the origin, which is a 'coincidence' in the 
CEM. 

The scenario discussed above is theoretically consistent 
with the Color Octet Mechanism (COM) |23|| , which is 
based on an NRQCD expansion in powers of the 
bound state velocity v. We only discussed the lowest 
order contribution in v, but added rescattering effects 
which there appears to be no reason to exclude in a COM 
approach. From a phenomenological point of view the 
scenario presented in [gj is quite different from the COM, 
in that the numerically important effects are ascribed to 
hard rescattering rather than to relativistic bound state 
corrections. 

A simple and explicit model for quarkonium produc- 
tion according to our scenario gave some encouraging re- 
sults that were not anticipated, in particular concerning 
the (non-)polarization of the J/ip and the cr(Xci)/c(Xc2) 
ratio. However, the model parameter p, which describes 
the spatial expansion of the heavy quark pair between its 
creation and rescattering times, had to be rather large 
(p ~ 3) to explain the measured a{Xcj)/&{J/ip) ratio. 

The hard rescattering scenario can also be applied to 
quarkonium production at large pj_, where the subpro- 
cess is gluon fragmentation to heavy quarks. It would be 
useful to make a quantitative study of this process in the 
framework of an explicit model. 

The scenario discussed above does not address the spe- 
cial features of quarkonium production observed at large 
Xf (increased nuclear suppression and longitudinal J/ip 
polarization). Other effects, e.g., as discussed in Ref. [ p5[ , 
need to be taken into account in that kinematic region. 
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